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M.  HELLER  and  J.  PAUL 


SUMMARY 

At  Materials  Research  Laboratories  (MRL)  Melbourne,  a  compact 
specimen  design  has  been  developed  which  is  suitable  for  both  plane- 
strain  fracture  toughness  and  Jjq  testing.  Compliance  and  stress 
intensity  factor  calibrations  are  given  for  the  new  design  of 
specimen  using  finite  element  analyses  and  the  results  are  compared 
with  those  for  the  ASTM  compact  specimen. 
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NOTATION 


Length  of  crack  measured  from  load  line 

Thickness  of  specimen 

Normalised  compliance 

Young's  modulus 

Critical  value  of  J  integral 

Stress  intensity  factor 

Plane  strain  fracture  toughness 

Normalised  stress  intensity  factor 

Load  applied  to  specimen  at  pin  hole 

Distance  from  crack  tip 

Radius  of  pin  hole 

Displacement  in  y  direction  of  node  on  crack  face  at  distance  r 
behind  crack  tip 

Poisson's  ratio 

Width  of  specimen  measured  from  load  line 
Cartesian  co-ordinate  axes  system 


(1) 
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1 .  INTROOUCTION 


The  fracture  mechanics  parameters  Kjc  Uthe  plane  strain  fracture 
toughness)  and  (the  critical  value  of  the  jj  integral)  characterise  the 
resistance  of  a  material  to  cracking  for  plastic  and  elastic-plastic 
conditions  at  the  crack  tip  respectively.  Standard  specimen 
configurations  are  available  for  and  Jj-Q  testing,  ajjd  these  are  given 

in  ASTM  standards  E-399  and  E-813  respectively  [  1 , 2] .  At  Materials 

- 

Research  Laboratories  (MRL)  Melbourne  a  compact  specimen  design  has  been 
developed  which  has  been  found  to  be  suitable  for  both  and 

testing.  For  determining  a  stress  intensity  factor  (K)  calibration 

is  required  for  the  specimen,  and  for  determining  Jjq  using  the  single 
specimen  technique,  a  compliance  calibration  is  needed.  This  memorandum 
describes  these  calibrations,  which  were  obtained  using  finite  element 
analyses,  and  gives  equations  which  have  been  fitted  to  the  results  for 


convenience. 
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2.  ANALYTICAL  DETAILS 

2.1  Specimens  and  Analysis 

The  ASTM  standard  specimen  configurations  used  for  Kjq  testing  are 
shown  in  Figs.  1  and  2.  Two  variants  of  the  MRL  compact  specimen  design 
were  investigated  and  are  shown  in  Figs.  3  and  ^4.  The  two  specimens  are 
identical  except  that  one  has  a  longer  notch  than  the  other.  The  short- 
notch  specimen  is  designated  as  specimen  A,  and  the  long-notch  specimen  is 
designated  as  specimen  B  (the  MRL  designations  are  J03  and  J*l3 
respectively).  Finite  element  methods  were  used  to  determine  specimen 
deformations  for  two-dimensional  conditions.  The  analyses  were  carried 
out  for  plane  stress  although  plane  strain  specimen  compliances  can 
readily  be  determined  from  the  plane  stress  values. 


2.2  Crack  Tip  Modelling 

Accurate  calculations  of  normalised  compliance  (C)  and  normalised 
stress  intensity  factor  (K)  require  detailed  modelling  of  crack  tip 
behaviour.  A  number  of  finite  element  methods  can  be  used  to  model  two- 


(2) 


dimensional  crack-tip  behaviour  and  several  of  the  better  ones  are 
reviewed  in  [3]  and  [A].  In  the  present  analyses  small  rectangular  iso 
parametric  elements  at  the  crack  tip  were  used  with  their  mid-side  nodes 

y 

shifted  to  the  quarter  points  [5]-  These  elements  give  the  required  r* 
displacement  singularity,  they  are  accurate,  and  are  easy  to  implement. 


2.3  Pin  Loading 

Both  specimen  geometies  were  analysed  for  both  point  and  distributed 
pin  loading.  For  point  loading  the  load,  P,  was  applied  to  one  node  at 
the  top  of  the  hole,  and  for  distributed  loading  the  load  was  spread 
evenly  across  nodes  at  the  top  of  the  hole  over  a  distance  equal  to  the 
hole  radius,  R. 


2.4  Noraalised  Conpllance 

Values  of  normalised  compliance  were  determined  using  the  equation, 
given  in  [2],  namely 


2EBV 

P 


(1) 


where  V  is  the  displacement  of  the  sharp  corner  on  the  specimens  in 
the  loading  direction  (see  point  marked  x  in  Figs.  3  find  4.) 


2.5  Nornalised  Stress  Intensity  Factor 

To  determine  stress  intensity  factors  for  the  various  specimen 
geometries  the  equation  given  in  [6]  was  used,  namely 


K  - 


uE  2n 

4(1-v^)  /  ’’ 


(2) 


(3) 


where  u  is  the  y  displacement  of  a  near-tip  node  on  the  crack  face 
behind  the  crack  tip  and  r  is  the  distance  of  that  node  from  the  crack 
tip.  Since  equation  (2)  is  valid  only  very  close  to  the  crack  tip,  a 

corner  node  (eg  node  '2')  of  the  crack  tip  element  was  used  in  the 

analyses  as  indicates  schematically  in  Fig.  5. 

Normalised  stress  intensity  factors,  K  ,  were  then  obtained  using 
the  equation 

K  =  -  ■  (3) 

where  w  is  the  specimen  width  measured  from  the  load  line.  This 

equation  is  used  since,  for  a  compact  specimen  at  a  given  crack  length 

a  P 

ratio,  -  ,  K  is  proportional  to  -  [?]. 

B/w 

3.  NUHERICAL  ANALYSES 

Finite  element  analyses  were  made  with  the  PAPEC  suite  of  programs  on 
the  ARL  VAX  11/780  computer.  The  stiffness  matrices  were  computed  using  2 
X  2  reduced  integration  and  all  solution  steps  were  performed  using  double 
precision . 

3.1  Specimen  A 

The  elastic  properties  of  the  material  were  taken  as  E  =  210  GPa  and 
V  =  0.3.  The  specimen  thickness  was  B  =  25.5mm  and  a  tensile  load  of  P 
=  4000N  was  applied  at  the  pin  hole. 

Analyses  were  conducted  for  values  of  crack  ratio  -  between  O.36 

and  0.76,  for  both  point  and  distributed  loading.  A  typical  finite 
element  mesh  with  point  loading  is  shown  in  Fig.  6(a) 

(for  =  0.52)  and  the  same  mesh  with  distributed  loading  is  shown  in 

Fig.  6(b).  The  mesh  consists  of  15^  eight-noded  isoparametric 

quadrilateral  elements  and  9  six-noded  isoparametric  triangular 
elements.  As  discussed  in  Section  2.1  the  crack  tip  elements  have  their 


mid-side  nodes  shifted  to  quarter  point  positions,  and  these  elements  were 
typically  th  the  length  of  the  crack. 


3.2  Specimen  B 

The  elastic  properties  of  the  material  were  taken  as  E  =  73  GPa  and 
V  =  0.3.  The  specimen  thickness  was  B  =  6mm  and  a  tensile  load  of  P  = 
632ON  was  applied  at  the  pin  hole. 

Analyses  were  conducted  for  values  of  crack  ratio  -  between  0.A6 

w 

and  0.7,  for  both  point  and  distributed  loading. 

A  typical  finite  element  mesh  with  point  loading  is  shown  in  Fig. 

7(a)  (for  =  0.52),  and  the  same  mesh  with  distributed  loading  is 

w 

shown  in  Fig.  7(b).  For  this  specimen  the  mesh  consisted  of  212  eight- 
noded  isoparametric  quadrilateral  elements  and  11  six-noded  isoparametric 
triangular  elements. 

Again  the  crack  tip  elements  had  their  mid-side  nodes  shifted  to 
quarter  point  positions,  and  for  this  specimen  these  elements  were 
typically  — ^  th  the  length  of  the  crack.  (The  mesh  for  this  specimen 
was  more  highly  refined  than  strictly  necessary  because  it  was  being  used 
for  other  work  involving  plasticity  analysis). 


i(.  RESULTS 

Values  of  normalised  compliance  (C)  calculated  using  equation  (1) 
for  specimens  A  and  B  are  given  in  Tables  1  and  2.  For  any  value  of 
considered  the  difference  in  C  values  between  specimens  or  loading 
is  very  small  and  may  be  neglected  when  used  for  test  purposes,  for 
example,  as  can  be  seen  by  comparing  values  at  =  0.56  or  0.60.  A 

least-squares  best  fit  polynomial  equation  was  determined  for  the  poolled 
data  over  the  range  0.il6  <  <  0.7,  with  the  following  result. 


C  «  b  + 


b  (-)  +  b 
1  u  2 


,a  * 

(-)  ♦  b 

w  * 
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(4) 


(5) 


where 


“3630.76 

33917.16 

“125878.0 

233^6.15 

-216373.5 

80m2.6^ 


Equation  (4)  is  within  +  0.35?  of  the  finite  element  results  for  any 
point  within  the  range  noted  above. 

It  should  be  noted  that  values  of  normalised  compliance  for  plane 
strain  conditions  can  be  determined  readily  from  the  present  results  using 
the  equation  given  in  [7]: 


C  =  C  ( 1 

plane  strain 


v^) 


(5) 


Values  of  normalised  stress  intensity  factor  (K)  calculated  using 
equation  (3)  for  specimens  A  and  B  are  given  in  Tables  3  and 

3  — 

For  any  value  of  -  considered  the  difference  in  K  values  between 

w 

specimens  or  loading  is  also  very  small  and  may  be  neglected  when  used  for 

test  purposes  as,  for  example,  can  be  seen  by  comparing  values  at  -  = 

w 

0.56  or  0.60.  Again,  a  least-squares  best-fit-polynomial  equation  was 

fitted  to  the  poolled  results  over  the  range  0.46  <  -  <  0.7.  The 

w 

derived  equation  is 
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(2  + 


(1 


a  ,3/2 
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»2  * 


(6) 


where  b^ 


b 

b 


2 
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0.21291 

6.682871 

-12.29665 

7.320672 


Equation  (6)  is  within  0.35J  of  the  finite  element  results. 

Averaged  values  of  C  and  K  (determined  from  finite  element  analyses) 
for  specimens  A  and  B  are  compared  in  Tables  5  and  6  with  those  for  the 
ASTM  compact  specimen  for  selected  values  of  The  C  values  for  the 

ASTM  specimen  were  obtained  from  [2],  and  the  K  values  were  evaluated 
using  equation  (7.^)  given  in  [7]  in  conjuction  with  equation  (3)  of  this 
memorandum.  It  can  be  seen  that  the  MRL  specimens  have  slightly  higher 
values  for  both  C  and  K  than  the  ASTM  specimens. 

Since  the  values  of  C  and  K  determined  for  both  the  ASTM  and  MRL 
specimens  and  given  in  Tables  5  and  6,  would  be  accurate  to  better  than 
0.5X  the  differences  in  C  and  K  shown  in  both  Tables  5  and  6  respectively 
can  be  attributed  to  specimen  configuration  and  not  to  calibration 
procedure. 

5.  CONCLUSION 

Normalised  compliance  (C)  and  normalised  stress  intensity  factor 
(K)  calibrations  for  MRL  compact  specimens  A  and  B  have  been  carried  out 
using  finite  element  analyses. 

Polynomial  expressions  have  been  fitted  to  the  C  and  K  calibration 

results  over  the  range  0.^6  <  -  <  0.7. 

w 

The  two  specimens  have  essentially  the  same  values  for  the  parameters 
C  and  K  over  the  range  of  crack  lengths  considerea  even  though  they  have 
different  notch  lengths.  The  MRL  specimens  however,  have  higher  values 
for  both  C  and  K  than  the  ASTM  compact  specimens  due  to  differences  in 
geometry. 
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Specimen  A  for  Plane  Stress  Conditions 


Normalised  compliance 


Point  pin  loading 


Distributed  pin  loading 


TABLE  2 

Normalised  Compliance  Value  for  MRL 
Specimen  B  for  Plane  Stress  Conditions 


a/w  Normalised  compliance 


Point  pin  loading 

Distributed  pin  loading 

0.46 

34.49 

34.63 

0.48529* 

38.46 

38.60 

0.50882* 

42.35 

42.92 

0.56 

54.68 

54.86 

0.60 

67.59 

67.72 

0.66 

96.71 

96.82 

0.70 

127.43 

127.51 

used  for  additional  plasticity  analysis 


a/w 

Normalised  stress 

intensity  factor  K 

Point  pin  loading 

Distributed  pin  loading 

0.i(6 

1 .955 

1 .954 

0.^8529 

2.128 

2.127 

0.50882 

2.293 

2.292 

0.56 

2.7*11 

2.740 

0.60 

3.195 

3.194 

0.66 

4.160 

4.158 

0.70 

5.106 

5.106 

»  gA  A.*  i.»  / 


MRL  Compact  Speciaena  and  ASTM  Compact  Specimen 


Normalised 

Compliance  C 

Percentage 

difference 

MRL  compact 

ASTM  compact 

specimens  * 

specimen 

0.50882 

.56 

.60 

.66 

.70 


^*2.35 

5‘t.68 

67.59 

96.71 

127.43 


39.71 

50.7 

63.5' 

92.1 

122.5 


»  Finite  element  results  for  specimens  A  and  B  averaged. 

TABLE  6 

Comparison  of  Normalised  Stress  Intensity  Factor  Values  for 


Normalised  stress^ 

intensity 

Percentage 

factor  K 

difference 

MRL  compact 

ASTM  compact 

specimens  * 

specimen 

0.50882 


Finite  element  results  for  specimens  A  and  B  averaged. 
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On  page  3  the  equation. 
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is  incorrect.  It  should  be. 
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